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Abstract 

A macroscopic elastic description of stresses in static, preloaded gran- 
ular media is derived systematically from the microscopic elasticity of in- 
dividual inter-grain contacts. The assumed preloaded state and friction at 
contacts ensure that the network of inter-grain contacts is not altered by 
small perturbations. The texture of this network, set by the preparation of 
the system, is encoded in second and fourth order fabric tensors. A small 
perturbation generates both normal and tangential inter-grain forces, the 
latter causing grains to reorient. This reorientation response and the in- 
cremental stress are expressed in terms of the macroscopic strain. 

The transmission of stress in granular media has a rich phenomenology 
D 121 nil as illustrated by the emblematic sand pile problem. In a conical pile 
obtained by pouring grains from a point source (hopper outlet), the pressure 
profile at the base of the pile is not proportional to the height of the pile, as it 
would if the weight of each grain were transmitted strictly vertically; nor does 
the pressure vary monotonically from edge to center, as predicted by traditional 
isotropic elastic or elasto-plastic approaches. Rather, the pressure profile devel- 
ops a local minimum (often termed 'pressure dip') at the center of the pile base, 
below the apex of the pile 01|Slini[7]. By contrast, if a pile of the same shape is 
prepared layer by layer, e.g., by sprinkling from an extended source, the pres- 
sure profile does acquire a maximum at the center of the base ('pressure hump') 
m m C] . These and similar experiments indicate that the local structure 
of the pile (often called 'texture'), which governs stress transmission, depends 
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sensitively on the preparation of the system. In two dimensional packings of 
monodisperse disks |7], for example, the observed distribution of inter-grain 
contact orientations reflects the local grain ordering, with a clear six-fold mod- 
ulation. However, as expected from simple symmetry arguments (see Fig. ^ 
for an illustration), the distribution displays a stronger degree of anisotropy 
(including two- and four-fold modulations) in a usual pile, built with a point 
source, than in a pile grown by sprinkling grains. Further experiments [3|H|, as 
well as simulations |5] , confirm that the anisotropy of the microscopic contact 
distribution is correlated with the presence of a macroscopic pressure dip. 

As it proved difficult to fit the wealth of observed stress proffies IHI^ISllTIITni 
imiS) with usual elastic theories, in the past decade physicists introduced highly 
ideahzed discrete models of static granular media 0)^1 El El based on proba- 
bilistic rules of force transmission between neighbouring grains. Although these 
models reproduce experimental features, they leave stresses underspecified (and 
consequently rely on a somewhat ad hoc constraint among stress components), 
as a result of an incomplete treatment of inter-grain forces. Furthermore, care- 
ful investigations of force transmission models |15l I16L ITT] seem to indicate that 
disorder in the granular packing might cause it to behave again as an elastic 
medium at large scales. In this spirit, Ref. El reviews the various possible out- 
comes of an (orthotropic) El anisotropic elastic theory [201, without however 
addressing the origin of the anisotropy in the context of granular media. 

In a continuum, any volume element is compactly embedded in the surround- 
ing medium. By contrast, in more fragile granular media, stress is transmitted 
through isolated inter-grain contacts, and under a perturbation grains slightly 
reorient with respect to their surroundings. This occurs even in the case of 
a simple compressive strain: consider for example a compression along the x- 
axis and focus on the contact of a given grain with a neighbouring grain. The 
torque that this contact imparts to the grain vanishes only if the contact is 
parallel or perpendicular to the x-axis. Thus, generically, a non-vanishing total 
torque is applied on a grain by contacting grains, and consequently it reorients 
so as to restore the torque to zero in equilibrium. A continuum theory derived 
from inter-grain forces must include effects of grain reorientation. In this vein, 
a generalization of usual elasticity that accommodates the 'micro-rotation' of 
points in addition to the 'macro-rotation' of the medium was formulated by the 
Cosserat brothers and subsequently generalized and applied to 'complex' 
continuum media [^1^ . 

In the present work, our starting point is the well-established elastic theory 
of inter-grain contacts |24[I25| : on this basis, we construct a specific macroscopic 
elastic theory El that encodes the texture of the granular network. We consider 
a preloaded granular medium |27| and small applied incremental stresses, so that 
the corresponding elastic response be linear. While the unperturbed inter-grain 
forces reflect the rolling and sliding El that occur during the preparation of 
the system, as well as the external preloading, if any, we expect the response 
of the medium to small incremental stresses to be sensitive only to properties 
of the (unperturbed) contact network. Below, we derive general expressions for 
the (anisotropic) elastic linear response and the associated grain reorientation 
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Figure 1: In a pile constructed layer by layer (top), the local texture of the 
packing has only one prefered (vertical) direction and is therefore statistically 
symmetric about any vertical plane (and hence invariant under rotation about a 
vertical axis). By contrast, in a usual pile grown with a point source (bottom), 
the presence of another special direction (the downhill direction) lowers the 
symmetry, and the local texture is statistically symmetric with respect to a 
single vertical plane, that which contains the downhill direction. 

field. 

The contacting region between two preloaded spherical grains A and B (Fig.J^a)) 
is a disk with diameter a « [3FD{1 — i''^)/8E]-^/^ , where F is the inter-grain 
compressive force, D the grain diameter (or, equivalently, the centre-to-centre 
distance) , E the Young modulus and the Poisson ratio of the grain consti- 
tutive material . Because of the high friction established by the preloading, 
essentially no slip occurs at the contact upon a perturbation of the grains |29| . 
Consequently, the linear elastic response of such a contact is characterized by 
stiffness constants for compression or traction along the centre-to-centre vec- 
tor D (henceforth called contact vector) (Fig. J^b)) and k± for transverse forces 




Figure 2: Schematic illustration of the preloading and deformation of two con- 
tacting grains, (a) The diameter a of the contact region between two preloaded 
grains is non- vanishing (over-emphasized on the figure). A restoring force about 
this compressed state is generated by (b) further compression (or traction), (c) 
shear, or (d) reorientation. If a <^ Z), (c) and (d) are equivalent for small 
perturbations. Correspondingly, tensors e and fl enter Eq. (0) on the same 
footing. 
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due to shearing (Fig.J^c)) Transverse forces are also generated by the re- 
orientation of two contacting grains (Fig. J^d)). Typically, a <^ D, implying 
fc||, k± ~ Ea, and the ratio k±/k\\ is fixed by for example, for a purely com- 
pressively preloaded contact k±/k^^ = (2 — 2z^g)/(2 — i/g) |2S]. In addition to 
compressive and shearing forces, torques may be transmitted through twisted 
inter-grain contacts, with a corresponding stiffness kt ~ Ea^ relating the twist 
angle to the torque. As long as a <C -D, such twisting torques contribute negli- 
gibly to the behaviour of the medium and we neglect them in the present work. 
(In the case of non-spherical ( e.g., faceted) grains, in which twisting modes may 
be more relevant, one can keep track of them by defining a torque flux tensor, 
often termed 'couple stress' in Cosserat theories |^E2, and proceed as with 
the stress tensor defined below.) Combining these various elements, we write 
the force exerted by grain B on grain A as 

/^^ - kff {uf - uf) - kf^^'^^Hl^Dt^, (1) 

where 

and the sum over repeated indices is understood. The tensor fl^ is antisym- 
metric and is defined such that the rotation matrix [exp (fi*^)]^^. « 6ij + 17^ 
describes the reorientation of grain G. We note that the mean reorientation, 
{n^ + n^)/2, results in a sheared contact (Fig. J^d)) and thus contributes to the 
force (Eq. (^), while the difference Sl^ — corresponds to grain rolling, which 
does not affect the force between spherical grains . 

In order to build a continuum description of a preloaded granular medium, 
we identify the displacements if' of the centers of the grains (labelled by G) with 
a smoothly varying, macroscopic displacement field u{f). Similarly, we identify 
the grain reorientation il'^ with a smooth reorientation field fl{r), to be deter- 
mined. In doing so, we neglect the non-affine component of the grain displace- 
ments. This approximation is exact in the limit of a simple crystalline granular 
network, in which all grains are equivalent, but we expect it to be good even away 
from this limit. If the packing disorder becomes too important, the medium can 
still be described by an elastic theory, although its effective stiffness constants 
may depart significantly from those obtained below ^Ij. Grain displacements 
can now be expressed in terms of the strain tensor, as uf — uf" = EijD^^ , where 
Sjk = {djUk + dkUj) /2 as usual, and the force transmitted through a contact 
whose contact vector lies along the direction a reads 

Ma) = hj (a) [Sjk - ^jk) Dk{a), (3) 

where kij[a) and Dk{a) refer to a contact normal to the direction a. We 
emphasize that, as it appears explicitly in Eq. |(2Jl, the grain reorientation Q. 
affects inter-grain forces; it does not represent the antisymmetric part of the 
displacement gradient (which corresponds to solid body rotation and does not 
generate stress). 
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For a description in terms of stresses, it is convenient to define the average 
number /^(a) da of contact vectors that he within a sohd angle da about di- 
rection a, per unit volume, and to introduce local fabric tensors |21 Q and P 
whose slow spatial variations reflect those of D, kj_, and /i, as 

Qkj ^ J k^DkDj ^ da (4) 

and 

Ptjki = / (fc|| - fc±) Ai da. (5) 

In equilibrium, the total torque J2g' Z^*^ ^ D'^'^ imparted to grain G by its 
neighbours (labelled by G') vanishes. Locally averaged in space, this condition 
becomes 



/■ 



/(a) X D{a)ii{a)da = (6) 
or, equivalently, 

QQ + qn^sQ- Qe. (7) 

This identity will be useful in the derivation of the stress in terms of the strain 
only, and confirms that grains reorient unless their contacts lie along a principal 
axis of the strain (on average), i.e., unless e and Q commute. In order to 
compute the stress, we note that (a)/i(a)da is the number of contact vectors 
per unit surface area, oriented within da about direction a, that intersect a 
surface normal to direction j. Since fi{a) is the j-th component of the force 
transmitted along such contact vectors, the stress tensor can be written as 

fi{a) Dj{a)^.{a)Aa. (8) 

While this expression is not symmetric in general, it is symmetrized as expected 
by the above vanishing-torque condition (Eq. {Tj)). Combining Eqs. (OIHI), we 
obtain 

/•OO 

n= e-'Q{eQ -Qe)e-'^ds or fl.j = ~ e^j, (9) 



and 



Qi 



e—^ 2QeQ e~'^ds + P : e or = 2-^^ e,, + P.jueki- 



lo 1^ + Ij 

(10) 

The right-hand expressions refer to a basis in which Q is diagonal, and qi is the 
i-th eigenvalue of Q; indices i and j are not summed over. As expected, the 
reorientation, 17, depends on fcj^ (through Q) but not on fc||. Equation ^ is 
qualitatively corroborated by the experimental observation (32) that f2 is largest 
in regions of large deformation. 
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Equations H9I1U|) constitute our central result and describe the response of 
a preloaded granular medium to an incremental strain. If the medium is statis- 
tically isotropic (/i = p/A-K constant), we recover standard isotropic elasticity, 
with Lame coefficients Al = pD'^{k\^ — A:^)/15 and /xl = pD'^{k\\/lb + k^/lQ), in 
agreement with the results of Chang and Gao |23[I26| . In particular, if the Pois- 
son ratio fg of the grain constitutive material is positive, that of the preloaded 
granular medium, v = Vg/(\Q — 6i/g), lies between and 1/14. (We recall, how- 
ever, that in the isotropic limit our coarse-grained theory is only approximate). 
In an anisotropic medium, the fabric tensors Q and P encode the texture of the 
medium through the dependence of D, k\\, k±, and fi on a. Indeed, if jj, varies 
with a or if the preloading stress is anisotropic, stiffness constants may also vary 
with a as contacts are then more or less compressed depending on their orienta- 
tions and may be sheared. The elastic response described by Eq. ((10)1 involves 
only the first few multipolar components of the quantities D^k±fi (through Q) 
and -D^(A:|| — k±)fi (through P). Monopoles contribute to the isotropic part of 
the response. Dipoles are absent due to the fore-aft symmetry of the contact 
distribution. Quadrupoles in Q and P and octupoles in P contribute to the 
anisotropic elastic response. Higher order multipoles play no role in the present 
linear theory. As to reorientation, given by Eq. ©, the response is affected 
solely by the quadrupolar component of D'^k^p,. 

Before concluding, we mention that it is possible to construct an iterative 
scheme for calculating stresses in a sand pile. It relies on the assumption that 
the mechanical noise due to avalanches of grains at the free surface, and other 
factors such as temperature fluctuations, do not cause signiflcant rearrangements 
of the contact network [3^. A rigorous implementation of an iterative scheme 
still requires knowledge of the complicated mechanisms involved in the growth of 
the pile, which include both grain avalanching [33] and the non-trivial frictional 
mechanics putting grains to rest |28[ I34| . Indeed, these determine both the 
angular contact distribution and the tangential compressive stresses at the free 
surface. 

As a first attempt at a quantitative understanding, one may focus on the 
simpler problem of calculating the macroscopic response to an applied point 
force j^Eni^SEllSSlEI for different possible textures of the granular medium. 
Using Eqs. 10 llOfl . we find that stress transmission is very sensitive to the 
texture as summarized by tensors Q and P, and comes in qualitatively diverse 
forms. Specifically, depending on the degree of anisotropy, the incremental stress 
profile may be single- or double-peaked jSS]) in agreement with Refs. [HI I18| : 
furthermore, the width of the peaks depends also on the degree of anisotropy 
and, in the single-peak case, may depart from its isotropic counterpart |35| . 

In sum, we have derived a macroscopic formulation of stresses in granular 
media. This formulation differs from earlier ones in that it incorporates the 
possibly anisotropic texture of the granular network as well as the reorientation 
of grains induced by macroscopic deformations. Moreover, the central objects 
in our formulation, the fabric tensors, are defined in terms of the microscopic 
parameters that characterize inter-grain contacts. 

One can wonder whether this description might bear some signature of a 
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salient experimental fact: the existence of highly stressed regions known as 'force 
chains' pi llOl lS). arranged in a percolating, filamentous network that appears to 
convey a large fraction of the stress through the medium. Analogously [3^1 EH 
to force chains, stresses in a crumpled elastic sheet are confined mostly within 
narrow regions, the folds, if the system is allowed to bend in the third dimension 
(of the embedding space) In our case, grain reorientation may be viewed 

as an additional degree of freedom which the system uses to relieve stress. It 
would be interesting to investigate whether this reorientational freedom of grains 
may favor, in a heterogeneous system, stress condensation reminiscent of force 
chains. 
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